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LHC has found hints for a Higgs particle of 125 GeV. We investigate the possibility that such a 
particle is a mixture of scalar and pseudoscalar states. For definiteness, we concentrate on a two 
Higgs doublet model with explicit CP violation and soft Z2 violation. Including all Higgs production 
mechanisms, we determine the current constraints obtained by comparing h — >■ 77 with h — > VV* , 
and comment on the information which can be gained by measurements of h — ¥ bb. We find bounds 
S2I < 0.83 at one sigma, where | S2 1 = (| S2 1 = 1) corresponds to a pure scalar (pure pseudoscalar) 
state. 

PACS numbers: 12.60.Fr, 14.80.Ec, 14.80.-j 



I. INTRODUCTION 

Recently, ATLAS [T] and CMS have reported tantalizing signals for a 125 Higgs particle. The exact properties 
of this particle will be probed in the coming years. In particular, it could be a mixture of scalar and pseudoscalar 
states. Some authors have already studied the possibility that the 125 GeV particle is a pure pseudoscalar state, both 
within the two Higgs doublet model (2HDM) [21 Sj and in a more general context [5]. However, as the authors point 
out, this possibility is at odds with the current h — > VV* signal, since a pure pseudoscalar state does not couple to 
VV (where V = Z. W). Thus, the bounds on h — > VV* should allow us to constrain the amount of the pseudoscalar 
component in the 125 GeV Higgs. In this article, we study this issue in the context of a 2HDM [5] with explicit CP 
violation and soft-breaking of the usual Z 2 symmetry. This model has been advocated in Refs. [THT4"]. 

We concentrate on models of type I, where all fermions couple to the same Higgs field, and models of type II, 
where the up type quarks couple to one Higgs field, while the down type quarks and the charged leptons couple to 
the other. The constraints placed by current data on the type I and type II CP conserving models (as well as on 
the lepton-specific and flipped models), have already been studied in Ref. [15], assuming that the 125 GeV particle is 
the lightest scalar, in Ref. [16], under the hypothesis that the 125 GeV particle is the heaviest scalar, and in Ref. [3], 
assuming a pure pseudoscalar. These studies assume Higgs production exclusively through gluon-gluon fusion and 
assume that the remaining scalar particles have very large masses. In our study of scalar-pseudoscalar mixing, we 
remove these restrictions, considering also inclusive production through vector boson fusion, associated production of 
a scalar and a vector boson and bb — >• H production. We also improve on Refs. [3l[l5l[l6] by allowing any scalar masses 
and mixings consistent with experiment and with the theoretical constraints from positivity, unitarity, perturbativity, 
and the oblique radiative corrections. 

In section |H] we describe succinctly the 2HDM with explicit CP violation and soft Z2 violation which we will use as 
a concrete example of scalar-pseudoscalar mixing, relegating to appendix [A] the formulae we have used. In section [TTT] 
we show our main results, and we conclude in section [TV] 
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II. A SPECIFIC MODEL FOR SCALAR-PSEUDOSCALAR MIXING 

As a specific example of scalar-pseudoscalar mixing, we will study a model with two Higgs doublets (fa, fa), with 
explicit CP violation, and with soft-violation of the Z 2 symmetry fa — > fa, fa — > — fa. As far as we know, this model 
was first written by Ginzburg, Krawczyk and Osland [7]. It was later studied in detail in Refs. [7HT5]. Here, we will 
follow the notation of Arhrib et. al 14J, which has a very clear presentation of this model. For ease of reference, we 
collect here some of its most important characteristics. 



The Higgs potential of this model is 



M "T-ll | , ,2 TO 22, , ,2 m 12 ,f - ( m i 2 )* A, 



Ai 

' 2 l ' r±l 2 



il 4 + ^I0 2 | 4 + Hfa\ 2 \fa\ 2 + A 4 {^\fa) 



+ y(^ 2 ) 2 + f {^Ifaf, (1) 

where hermiticity forces all couplings to be real, except m\ 2 and A5. If the latter are complex and arg(As) 7^ 2arg(m 2 2 ), 
then there is explicit CP violation in the Higgs potential. In addition, we also allow for explicit CP violation in the 
Yukawa terms, leading to CP violation through the CKM matrix, as is needed to account for the current data on CP 
violation in the K and B systems. The m\ 2 terms constitute a soft violation of the Z 2 symmetry, which does not 
affect the renormalizability of the theory (the renormalization group equations of the quartic terms do not depend 
on the quadratic terms). Naturally, one can change the phases of the fa and fa fields; any physical observable must 
be rephasing invariant |17j . An overall phase corresponds to a global hypercharge transformation, and has no effect 
on the lagrangian; it can be used to render the vacuum expectation value (vev) of fa real. Rephasing fa can now be 
used to render its vev also real, 

(fa) = Vl /V2, (fa) = V2/V2, (2) 
substantially simplifying the minimization conditions, which become 

rn?! = —Re (m? 2 ) — + Ai v\ + A345 v\ 

m 22 = — Rc(to 12 ) h A9 v 2 + A345 v, 

Im (T7j,f 2 ) = viv 2 Im (A 5 ) , (3) 

where A345 = A3 + A4 + Re(As). With our conventions, v — \Jv\ + v\ = (v^G^) -1 / 2 = 246 GeV. Thus, v\ and V2 
depend only on tan (3 — v 2 jv\ . 

We denote by "Z 2 basis", the basis where the Higgs potential has the form in Eq. (JlJ and the vevs are given by 
Eq. pi), and we parametrize the fields in this basis by 



(4) 



V2^ L ' ll A-±J J \ V2^ 

The charged fields are changed into the mass basis by 

G + = cos /3 4>l + sin f3 ip 2 , 

H+ = -SU1/30+ +cos(3(p+. (5) 
We apply the same transformation to the imaginary parts of the neutral fields, 

G° = cos/3xi + sin/3x2, 

i] 3 = -sin/?xi + cos /3 X2, (6) 

but not to the real parts of the neutral fields [THIEI]- As shown below, this is done in order to keep a clean definition 
for the angles ai leading the neutral fields from the Z 2 basis directly into their mass basis 1 . G + and G° are the 



1 Recall that /3 is the angle leading from the Z2 basis into the Higgs basis. In the completely Z2 symmetric case, a — /3 is the angle leading 
from the Higgs basis into the mass basis, and (thus) a is the angle leading directly from the Z2 basis into the mass basis. 
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would-be Goldstone bosons and H + is already the physical charged Higgs field, with mass m#± . Finally, one needs 
to diagonalize the (squared) mass matrix for the neutral fields M 2 , whose components are 

This is achieved through an orthogonal transformation 

= R f m I , (8) 

\ h 3 J \V3 / 

such that 

RM 2 R T = Ai&g (ml, m 2 2l ml) , (9) 
and mi < 7772 < 777.3 are the masses of the neutral Higgs particles. The matrix R may be parametrized by [9] 

C\C2 S\C 2 s 2 

R= I -(cis 2 s 3 + S1C3) C1C3-S1S2S3 c 2 s 3 I (10) 
-C1S2C3 + S1S3 -(cis 3 + SiS 2 C 3 ) C 2 C 3 

with s.j = sinaj and Cj = cosa^ (i = 1, 2, 3). Without loss of generality, the angles may be varied in the intervals [DJ 

- tt/2 < ai < tt/2, -7r/2 < a 2 < ir/2, < a 3 < vr/2. (11) 

The lightest neutral Higgs particle (the putative 125 GeV state) is determined by the first row of R. If \s 2 \ = 0, 
then 773 does not contribute to hi, which is a pure scalar. Notice that, in this case, there may be CP violation due to 
mixing in the h 2 , h 3 states [2D]. This possibility will affect CP violating observables but not the current known data. 
If \s 2 \ — 1, then only 773 contributes to hi, which is a pure pseudoscalar. In this case there is CP conservation. We 
conclude that 

\s 2 \ =0 =>■ hi is a pure scalar, (12) 
I S2 1 = 1 hi is a pure pseudoscalar, (13) 

and I S2 ] is a measure of the pseudoscalar content of the lightest Higgs scalar. 

Given Eqs. the scalar sector depends only on the eight parameters (3, Re(m 2 2 ), Im(mf 2 )) ^1,2,3,4, and Re(A5). 
As suggested in Ref. [10) . these can be traded for 7774, m 2 , m H ±, ai 2,3, /3, and R,e(m\ 2 ). In this approach, 7773 is a 
derived quantity, given by 

m 2 = raj ^13(^12 tan/3 - i? n ) + m| R 23 {R 22 tan^ - R 21 ) 

ma " R 33 (R3i -i? 32 tan/3) ' [ ' 

The implementation of the Z 2 symmetry in the fermion sector means that each fermion type (up quark, down 
quark, and charged leptons) can couple only to one of the original scalar doublets. As an example, let us consider 
three down type quarks coupling exclusively to <pi . We start from the Yukawa lagrangian 

- C Y = gzA^njj + h.c, (15) 

where q£ = (pl, n h) is a left-handed doublet of the gauge group, having 3 components in family space, while 77 r has 3 
down type quarks, each a singlet under the gauge group. Ti is the 3x3 matrix of Yukawa couplings, and h.c. stands 
for hermitian conjugation. After SSB, the lagrangian has a piece involving the neutral component of <f>i, which may 
be written as 



1 , Vi + iXi 

Vl 



M 

n R + h.c. = d L M d d R + d L — (7/1 + ixi) d R + h.c. (16) 

Vl 



To obtain the last expression we have rotated the fields njj and 77^ into the mass basis (d^ and d R ), diagonalizing the 
matrix T\Vi/\/2 to obtain the (diagonal) down quark mass matrix M d . Now, we need 

771 ±ixi = (Rki T ispRkz) hk ± icpG , 

r\ 2 ±ixi = (R k 2±ici3Rk3)h k ±ispG°, (17) 
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Type I 


Type II 


Lepton 
Specific 


Flipped 


Up 
Down 
Leptons 


S/3 p "p 

^ + icg^ 

"P P »P 

^+icg^- 

sp P sg 


sp p up 
cp p cp 
C P p C P 


— ICg — kA 

sp p up 

^+icg^ 

sp p up 

^-isg 1 ^ 

C P p C P 


Rl2 • -R13 

Kix-isgEm 

cp p cp 
sp p sp 



TABLE I: Couplings of the fermions to the lightest scalar, hi, presented, for each case, in the form a + ib. 



which were obtained by inverting Eqs. ^ and (J8|, and where a sum over k = 1,2,3 is implied. The couplings with 
the neutral scalar particles are obtained substituting Eq. (17) in Eq. (16) to get 



d L — (Rki ~ ispRm) h k d R + h.c. = d L — ( — - is^— 75 ] d R h k . (18) 
We are interested in the couplings of the fermions to the lightest scalar hi. Comparing with the notation of the 



effective lagrangian (All in the appendix, needed for our calculations, we obtain 



— — , b = —sr — — (down type quarks couple to <f>i) . (19) 
A similar analysis leads to Table [T] Notice that S2 = implies i?i3 = 0, in which case all b coefficients in Table [T] 



vanish, confirming that hi couples to the fermions as a pure scalar, and vindicating Eq. (12). Similarly, | S2 1 = 1 leads 
to i?n = = -R12, and all a coefficients in Table |T] vanish, implying that hi couples to the fermions with 175, being a 
pure pseudoscalar, as stated in Eq. ( p~3| ) . 

In Model I, all quarks and charged leptons couple to the doublet 4>2 , and the corresponding coupling factors can be 
read from the first column of Table |lj In Model II, the up quarks still couple to <f>2 but the remaining fermions couple 
to <fii - the respective couplings are shown in the second column of Table [I] The Lepton-Specific and Flipped models 
discussed in Refs. [T5l[TS] would be obtained by taking for the charged leptons the opposite choice taken for the down 
quarks. We have checked that our results reproduce those for Model II included in Ref. |14j . 

Expanding the covariant derivative terms of the neutral scalars, the triple interactions of hi with WW and ZZ 



may be written as in Eq. (A3), with 

C = cpRu + spR 12 . (20) 

As expected, when \s2\ = 1, -R11 = = R12 and there is no tree-level coupling of the pure pseudoscalar hi to a pair 
of gauge bosons. From the Higgs potential, we may get the triple vertex of hi with the charged Higgs bosons as in 



Eq. (|A2J), with 

- A = cp [s|Ai4 5 + c^A 3 ] Rn + s p [c|A 24 5 + s 2 ^] R i2 + spcp Im(A 5 ) R i3 , (21) 
where A145 = Ai — A4 — Re(As) and A245 = A 2 — A 4 — Re(As). In the pure pseudoscalar limit, only the last term of 



Eq. (21 ) survives, showing that a pure pseudoscalar can only couple to a pair of charged scalars if there is explicit CP 
violation in the scalar potential, through Im(As). 

For a given set of input parameters, the effective couplings (a, 6, C, and A) discussed in this section can be used 
on the equations in appendix |A"| in order to find the production rates for hi and its decay rates into all final states. 
For each final state /, we define the ratio 

R = °{PP -» hi) BR(/ti -> /) 
f a{pp -> K)sm BR(/i -> /) SM ' 1 ' 

where cr(pp hi) is an inclusive production rate obtained by summing over production mechanisms of hi, and 
BR(/i! — > f) is the branching ratio of the decay of hi into the final state /. 

The production mechanisms included in a(pp — > hi) include: the usual gluon-gluon fusion processes (which involve 
loops with both top and bottom quarks); vector boson fusion (VBF) processes; so-called associated production 
processes, with a W or Z boson in the final state, alongside hi; and the bb — > hi process as well. An important point 
needs to be explained at this juncture: since hi is a mix of scalar and pseudoscalar states, its production mechanism 
will involve, likewise, a mix of production rates pertaining to a scalar and a pseudoscalar particle. For instance, a pure 
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pseudoscalar particle has no VBF or associated production mechanisms; on the other hand, due to different fermion 
couplings, the gluon-gluon fusion cross section is different for production of a scalar or pseudoscalar states of the 
same mass. The exact formulae for cr(pp — > hi) can be found in Appendix [A] Using those formulae we were able to 
express the production factors in terms of next-to-leading order cross sections for each of the mechanisms considered, 
calculated using HIGLU [5T] for the gluon-gluon cross section, bb@nnlo [25] for the bb process and reference [33] (and 
references therein), for the remaining processes, all multiplied by the adequate factors pertaining to the mixing of 
scalar-pseudoscalar states. 

The branching ratios BR(/ii — > /) also need to be computed considering the mixed nature of the hi state, and the 
relevant formulae can be found in Appendix [A] Notice that, since the branching ratio is the partial width divided 
by the sum of all decay widths, a parameter choice that affects, for example, hi — > bb will have an impact on i? 77 
because it affects the overall decay width. 

We are particularly interested on what one can learn from the current LHC bounds on i? 77 , Rzz, Rww, an d 
R b l. These are combined with other known constraints, including recent results from the Tevatron [Ml US]. In this 
analysis, we utilize the explicit bounds obtained for the several Rf from the experimental data by Espinosa et. al 
[2"6] . Specifically, we use their Table III, updated after Moriond 2012. Similar summaries of the experimental bounds 
have been obtained in Refs. [27H29) . Since the current errors are large, we are not interested so much in the precise 
values but rather in the qualitative features hinted at by current bounds. Surprisingly, we can already exclude very 
large regions of the 2HDM parameter space, and place constraints on the pseudoscalar content of the lightest Higgs 
particle. 



III. ANALYSIS AND RESULTS 



We performed an extensive scan of the 8-dimensional parameter space of the models we are studying. Our fixed 
inputs are v = 246 GeV and mi = 125 GeV. We then took random values of: the angles a±, ai and a-i in their 
allowed intervals of variation, specified in Eq. (11); tan/3 between 1 and 30; the masses mi (above the value of mi) 
and m H ±, the latter with values above 90 GeV; and Re(m 2 2 ), taken in the range from — 10 6 to +10 6 GeV 2 . 

The charged Higgs couplings to the fermions have exactly the same form as the ones in a softly broken Z^ symmetric 
two-Higgs doublet model. Therefore the bounds derived for the charged sector still hold for the CP-violating scenario 
under study. The LEP Rf, constraint (from Z — > bb) excludes values of tan/3 < 1 even in the CP violating scenario. 
One should note that although the expressions for R b change in the CP-violating scenario, the main contributions for 
low tan /3 come from the charged Higgs diagrams [111 130] . Constraints from B-physics, and particularly those coming 
from b — > S7 [3TT436] . have excluded a charged Higgs boson mass below 300 GeV in models type II and Flipped, almost 
independently of tan/3. Charged Higgs bosons with masses as low as 100 GeV are instead still allowed in models 
Type I and Lepton-Specific [3"TH3"5] . Finally b — > S7 implies tan (3 > 1 for a charged Higgs mass below approximately 
600 GeV. 

Once a given set of parameters is chosen, a value for m\ is computed from Eq. (14). If that value is positive and 
larger than m 2 , we then verify that all theoretical constraints on the potential are satisfied. The quartic couplings 
A1...5 are computed using the formulae in Eq. (B.l) from Ref. [T3]. One then verifies whether these quartic couplings 
obey the conditions which ensure that the potential is bounded from below and that unitarity and perturbativity are 
satisfied - Eqs. (2.20) and (2.21) from Ref. [14], respectively. Finally, one verifies that the potential's parameters are 
such that the current constraints on the S, T and U oblique parameters are obeyed. To do so, we have used the 
formulae presented in Appendix D of Ref. [6 . If all of these constraints are obeyed, the set of parameters (a "point" 
in parameter space) is deemed satisfactory and one then proceeds to calculate the branching ratios and production 
factors of the hi state. 



A. Model I 



In Model I, as was explained earlier, all fermions couple to the same doublet, 2 by convention. The couplings 
between hi and the fermions are affected by the a and b parameters shown in eq. ( Al ) which, for this model and the 
tan/3 we chose, tend to be smaller than 1. We generated over 270000 points for this model. 

In Fig. [TJ we show how the i? 77 — Rzz plane is filled by Model I, for three selections of points: points for which 
one has an hi state which is essentially a scalar (i.e., for which the angle ct2 obeys |s 2 | < 0.1); an h\ state with 
0.45 < \s2\ < 0.55; and points for which hi is essentially a pseudoscalar (|s2| > 0.83). There are several salient 
features. First, we notice that even the current rather loose bounds already kill large regions of the model's parameter 
space. In particular, Rzz tends to be smaller than 1 in Model I, for any values of the parameters of the model. 
Second, a large pseudoscalar component (|s 2 | > 0.83) is already excluded, both by the ZZ bounds and by the 77 
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1.2 



FIG. 1: Scatter plot in the R 11 — Rzz plane for the type I 2HDM. The color coded points have the following correspondences: 
green (light grey) means S2I < 0.1 (hi is mostly scalar), blue (black) means 0.45 < |sa| < 0.55, and red (dark grey) means 
S2I > 0.83 (hi is mostly pseudoscalar) . The yellow and light blue (light grey and grey) bands shows the current ATLAS and 
CMS bounds, from Ref. [2B]. 



bounds. That the ZZ bounds constrain s 2 was expected and was the primary motivation for this work. Indeed, when 
|s 2 | = 1, Eq. ( 20 1 leads to C = 0, guaranteeing that there is no coupling of h\ to ZZ nor to WW. What is new is that, 
in this model, the known experimental bounds and theoretical constraints make |sa| ~ 1 inconsistent with i? 77 even 
if Rzz vanished. A third important point (not clearly visible in Fig. [T]) is that the 0.45 < | S2 1 < 0.55 (blue/black) 
region extends further into the i? 77 experimentally allowed region than the \s%\ < 0.1 (green/light grey) region 2 . This 
has the following implication: if the current ATLAS central values for i? 77 remain as the errors get smaller, then the 
model seems to prefer a mixture of scalar and pseudoscalar components in the lowest lying Higgs particle. Finally, 
we find that only for | S2 1 < 0.83 do we start to get points inside the band allowed by current experiments. 

We have mentioned that there are very few points above around i? 77 = 1 corresponding to the (almost) pure scalar 
solution |s 2 | < 0.1. There are quite a few more points in the large R 11 region corresponding to 0.45 < |s 2 | < 0.55. 
Fig. [I] shows that the former tend to concentrate around Rzz ~ 0.8. A better handle on the difference between 
|s 2 | < 0.1 and 0.45 < |s 2 | < 0.55 is provided by R bb , as shown in Fig. |2j where Model I's results are presented in the 



i? 77 — R b i plane. We see that the pure scalar solution with large i? 77 leads to R b i ~ 1, while the large values of R. 



correspond to R bb ~ 0.5 when 0.45 < |s 2 | < 0.55. Thus, as the errors get smaller, a comparison between i? 77 , Rzz, 
and R bb can be used to further constrain |s 2 | in Model I. 

Regarding the current bounds on R bb , Espinosa et al |26] provide an LHC interval of R bb < 3.3, which spans the 
entire range for this variable displayed in Fig. [2j In that reference there is also a bound stemming from Tevatron 
data, 1.3 < rTEVATron < 2 g. once translated into LHC bounds for Model I - that is, with the appropriate LHC 
production factors - this would give approximately R bb > 1.2 in Fig. [2j This would eliminate most of the available 
parameter space left by the bounds on i? 77 . Accurate measurements of R bb are thus of extreme importance for this 
model. 

We have also looked at hi — > W + W~ . Notice that, within our tree-level calculations, and even for a mixed hi 
state, Rww — Rzz, since we are dealing with ratios to SM quantities. As such, one can glean information about 
Rww from the plot in Fig. [T] After an earlier central value compatible with the SM, ATLAS presented at Moriond 



2 Please notice that there is some superposition of points in several regions, so that "underneath" the points shown in blue/black, for 
instance, there may well exist green/light grey points. This is common to all figures shown in this work. 
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FIG. 2: Scatter plot in the i? 77 — R b i plane for the type I 2HDM. The color codes are the same as the previous figure's. 

2012 results which are now slightly over one sigma smaller than the SM - compare the respective entries in Tables I 
and III of Ref. [25]. Their Moriond 2012 bounds impose R ww smaller than 0.8 for ATLAS and 1 for CMS. If this 
situation were to remain as errors get smaller, then the SM itself would be in trouble. As seen in Fig. [I] for hi -> ZZ 
in Model I, the values for hi — > W + W~ are smaller than one, enabling a possible accommodation of a low Rww 
signal. 

B. Model II 

We now turn to Model II. As explained, in this model the down-type quarks and charged leptons couple to the <pi 
doublet, and the up- type quarks couple to 4>2- The corresponding couplings can be read off Table |TJ and we see that 
the bottom quark couplings are proportional to l/cos/3. As such, one can expect, in this model, an enhancement 
of the production of hi via the gluon-gluon fusion 6-quark triangle, or the bb —¥ hi channel for large values of tan ft. 
Similarly, the branching ratio BR{h\ — > bb) will also reach higher values than in Model I for large tan/3. As was 
mentioned earlier, in Model II there are stringent constraints on the value of the charged Higgs mass - it needs to be 
larger than about 300 GeV - stemming from b — > S7 data. We included that cut in our scan of this model's parameter 
space, and again we generated over 270000 points for this model. 

In Fig. [3] we show our results for Model II in the i? 77 — Rzz plane, while Fig. [4] shows the results in the i? 77 — R bb 
plane. As expected, large values of are excluded by the ZZ bound. There are similarities and differences between 
the two models. The striking feature that i? 77 by itself constrains |sa j is common to models I and II, and in both 
cases I S2 1 is found to have to be smaller than about 0.83. 

The most noticeable difference between both models is that Model I keeps roughly i? 77 < 2, Rzz < 1, and 
R bb < 2.5 while, in Model II, values as large as i? 77 < 2.5, Rzz < 2.7, and R b i > 2.5 3 are allowed. Also, in 
Model II, the | s 2 1 < 0.1 (green/light grey) region extends further into the i? 77 experimentally allowed region than 
the 0.45 < I S2 1 < 0.55 (blue/black) region. This is the opposite of what we observe in Model I. Finally, we see from 
Fig. |4| that the pure scalar solution with large i? 77 tends to imply R bb < 1, while the large values of i? 77 correspond 
to R bb > 1 when 0.45 < | S2 j < 0.55. Thus, as better measurements are available, R bb might become instrumental 



3 Values as high as 10 for R b ^ were found, for extremely small values of R 11 , though they are not displayed in the plot of Fig.[4]for ease of 
presentation. These high values correspond to both the gauge-phobic limit of this model and the enhancement with large tan /? alluded 
to earlier. 
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FIG. 3: Scatter plot in the _R 77 — Rzz plane for the type II 2HDM. The color codes are the same as previous figures'. 
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FIG. 4: Scatter plot in the i? 77 — R b i plane for the type II 2HDM. The color codes are the same as previous figures'. 

in constraining s%. In fact, the current Tevatron bounds, cited by [33], translated into LHC for Model II, would 
give roughly 1.1 < R b i < 2.8, which, if confirmed, would exclude an extremely significative portion of this model's 
parameter space. In particular, it would seem to disfavor the possibility of hi being a pure scalar. 

Finally, a word on the constraints emerging from Rww for this model. As mentioned above for the Model I analysis, 
current ATLAS and CMS bounds favor values of Rww smaller than 1 - since Rww — Rzz, we can see from Fig. [3] 
that that region is heavily disfavoured for Model II by the current bounds on i? 77 . As a result, very low experimental 
values for Rww will exclude the SM and also the type II model. 
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IV. CONCLUSIONS 



With the LHC providing physicists with a wealth of data, and the first hints of the existence of a Higgs particle 
with a mass around 125 GeV, it becomes possible to constrain, not only the SM but also extensions of it, such as 
the 2HDM. Previous studies considered the possibility that the LHC might be observing a pure scalar, or a pure 
pseudoscalar. We have considered the possibility that the putative scalar candidate at the LHC is a mixed state, 
neither scalar nor pseudoscalar. We considered a specific model, a version of the 2HDM with a Z2 discrete symmetry 
which has been softly broken, such that the model has explicit CP violation in the scalar sector. The degree of 
"pseudoscalarity" of the lightest scalar h\ is measured by a mixing angle a.2 such that \si \ = | sinc^l — 1 corresponds 
to a pure pseudoscalar state (and \si\ — to a pure scalar one). We considered two specific extensions of the Z-i 
symmetry to the fermionic sector, the so-called Models I and II, which have very different phenomenologies. 

We have computed the production rate of hi times its branching ratio into several final states - Rf - relative to 
the expected values for such observables in the SM. We have shown that this version of the 2HDM can do at least 
as good a job as the SM in fitting the current data. But we have also shown that even the current loose bounds on 
i? 77 and Rzz already put severe constraints on the parameter space of these versions of the 2HDM. In particular, 
our work suggests that if current trends in the values of Rww a t the LHC, as well as in the Tevatron data on 
persist, the versions of the 2HDM herein considered might have a hard time reproducing the data. However, if that 
were the case, the SM would also be in trouble. Current data forces | S2 1 < 0.83, a constraint which holds roughly 
even if one considers only the i? 77 bound. This by itself excludes a large pseudoscalar component. 
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Appendix A: Production and decay Rates 



In this appendix we present the decay rates for a particle h with both scalar and pseudoscalar components. The 
relevant pieces of the Lagrangian are: 



C Y = - 



1 

-1 = - [VZGf+j 2 rrif ip (a + ib-f 5 ) ip h, 
£hH+H- = \vhH + H~, 

r+Tiru— 1 9 



c 



hVV 



= c 



gmwW+W* 1 - 



2cw 



h, 



(Al) 
(A2) 
(A3) 
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where a, b, and C are real, cw = cosOw, and 9w is the Weinberg angle. In the SM, a = C = 1, and b = A = 0. 
We find: 



T(h 77) = — - — 

V ; 128v/27T 3 



E 7V C Q ) a Ai /2 (r f ) + CA X (t w ) - ■ 7 ^ r - A A Q (r± ) 



J2N c Q 2 f bAf /2 (r f ) 
f 



(A4) 



— 1/2 

where v = [V^G^] « 246 GeV and mjj± is the mass of the charged Higgs. For ease of reference, we have used a 
notation close to that of Djouadi HQ] HI], where 



A 1/2 (r) = 2[t + (t-1)/(t)]i 
A? /2 (t) = 2r-V(r), 



and 



A)(r) 



/(r) = 



-[2r 2 + 3r + 3(2r- l)/(r)]r- 2 
-k-/(r)]r- 2 , 



[arcsin(y / r)] i 



T < 1 



T > 1 



The scaling variables are r, = Af 2 /(4m 2 ), where is the mass of the particle in the loop. 
Similarly, for the decays into two gluons we find 



(A5) 
(A6) 

(A7) 
(A8) 



(A9) 



(A10) 



where the sums run only over quarks q. 
The decays into fermions are given by 



r(/i -> //) = n c 



G/j. mj 
4v/27r 



M h [a 2 (3 3 f + b 2 (3 f ] 



(All) 



where /3/ = Jl — Arnj/Mj* = y/l — r 1 , while the decays into two vector bosons are given by 

v(h -> vWv w ) = c* 2 r SM (/i -> y ( *V ( * } ), 



(A12) 



and the partial decay widths in the SM-Higgs case in the two-, three- and four-body approximations, Tsy[(h —> 
V r Wy(*)) i can be found in Section 1.2.2 of Ref. HQ]. Notice that in no decay is there interference between the scalar 
a couplings and the pesudoscalar b couplings. 

The same is true in the production mechanisms. We find 



G a 2 

a(gg -> h) = 

512V27T 



(A13) 



where the sums, which run over all quarks q, are dominated by the triangle with top in the loop with, depending on 
tan /3, relevant contributions from the triangle with bottom in the loop. Therefore, we use 



a{gg^h) _ \a t A 1/2 {r t ) + a b A 1/2 (n)\ 2 + \b t Af /2 (r 4 ) + b b Af /2 (n)f 



a s ™(gg^h) 



\A 1/2 ( Tt ) + A 1/2 ( n )\< 



(A14) 
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Similarly, 



&VBF 
°VBF 



(J 



SM 
VH 



c 2 



and 



a(bb -> h) 
a SM (bb^h) 



(A15) 



(A16) 



Notice that these expressions hold for any model with the effective lagrangians of Eqs. ( Al |-( A3) 
Our results agree with those of Choi et. al \E 



% H3"]. For ease of reference, we show their definitions in terms of ours: 

F sf (r) = \A l/2 {r), 

F pf {r) = ±Af /2 (r), 

F 1 (t) = -Ax(r), 

F (r) = A (r), (A17) 



and 



y s f 



Ci — A, 



(c/30 2 ,i + apO^i) = C, 



-(\/2G M ) 1 / 2 m / a = 
-(VZG^y^mfb = - 



(A18) 
(A19) 
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